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Abstract 

The functional integral has many triumphs in elucidating quantum theory. 
But incorporating charge fractionalization into that formalism remains a challenge. 

This conference celebrates the achievements of path and functional integration in quantum 
physics. But it is good to remember that some of these achievements were hard to attain, 
because they required resolving unexpected subtleties of the functional formalism. Two 
historical examples will illustrate my point, and then I shall posit a new challenge. 

For the first example, we look at the description of rotationally symmetric motion in 
a potential V{r). The classical effective Hamiltonian reads H = + i^pi + V^(^"), where 
the centrifugal term represents r^^^/2m and vanishes for circularly symmetric motion. 
But in the quantum description of this motion, when it is confined to a plane, the classical 
centrifugal term is replaced by /i^(M^ — l/4)/2mr^, where M is any integer. In particular a 
residual attraction remains for s-waves, M = 0. This quantal attraction has the important 
physical consequence that in planar physics bound states exist, no matter how weakly 
attractive might be the potential V . (In this way, planar bound states follow the behavior 
of one dimensional bound states, rather than three dimensional, which require a minimum 
strength to achieve binding.) 

The functional integral involves integration over c-number functions. A direct change of 
variables from Cartesian to circular coordinates reproduces the classical centrifugal barrier, 
but misses the residual, attractive potential —f^j^mr^. This is particularly vexing, since 
coordinate changes are point canonical transformations, which are allowed in quantum 
mechanics. 

The resolution of this problem was given by Edwards and Gulyaev [Ij. But it required 
returning to the discretized formulation of the path integral and realizing that the angular 
step A0, is not of order of the temporal step Ar, but rather 0(V At). [This issue reappeared 
in the collective coordinate quantization of solitons. Investigations based on functional 
integrals missed 0(^^) terms that were found in an operator approach [2].] 
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For my second example, I turn to the anomaly phenomenon: the classical action can 
possess an invariance against a transformation that is not an invariance after quantiza- 
tion. Where is this effect in the functional integral, which involves "functional integration" 
over the exponentiated classical action? The answer, which was found by Fujikawa (after 
anomalies were discovered by conventional methods), located the effect in the functional 
measure. Evidently in the anomalous situation it is not invariant; a fact that is established 
after the measure is discretized [3]. 

In both instances we see that the validity of formal changes of variables in a functional 
integral must be assessed by a return to the discretized functional sum, in fact, by a return 
to conventional quantum theory. 

Now I shall present another peculiar quantum effect, which as far as I know has not had 
a functional integral description, even though a conventional quantal argument establishes 
it rather easily. I have in mind the phenomenon of fermion charge fractionalization in the 
presence of a topological defect [1]. 

It is well known and also guaranteed by various index theorems, that the Dirac equa- 
tion in the presence of a topological defect (kink in one spatial dimension, vortex in two, 
monopole in three etc.; we call these "solitons") possesses an isolated, normalized zero- 
energy, mid-gap bound state. Of course the Dirac equation also possesses positive energy 
solutions and negative energy solutions. The conventional instruction, given by Dirac, is 
to define the vacuum by filling the negative energy states and leaving the positive states 
empty; with this definition the vacuum charge vanishes. But what should one do with 
the mid-gap, zero-energy bound states, if it is present? Dirac is silent on this; he did 
not know about mid-gap bound states. The answer is that there is double degeneracy in 
energy, since filling the mid-gap state costs no energy. Moreover, the empty mid-gap state 
carries fermion charge —1/2, and the filled one, +1/2. Remarkably, this effect has been 
observed in polyacetylene — a 1-dimensional lineal material — and it has been proposed 
for 2-dimensional planar graphene. (Note: the Dirac equation is relevant to these con- 
densed matter materials not because of relativistic considerations; rather a well defined 
linearization of the energy dispersion near the Fermi surface gives rise to a linear matrix 
equation, which is of the Dirac type.) 

Fermion number fractionalization is established by the following argument. Consider 
the Dirac Hamiltonian H{ip) in a background field which can be topologically trivial 
or non-trivial. In the trivial case takes a homogenous value (fy, which gives rise to 
a mass gap in the Dirac spectrum, between the positive and negative continuum energy 
eigenstates, which we call "vacuum states" and denote them hy ipE- 



With the topologically non trivial background ip takes a soliton profile ips and the Dirac 
continuum eigenfunctions are called "soliton states," denoted by ^e- Also there is an 
isolated normalized state at zero energy, ^'q. 





H{ips)'^E 







(2) 
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We wish to compute the charge density in the presence of the sohton. This is defined 
relative to the charge density in the vacuum, where the background field is topologically 
trivial. 



pir) = J dE(^mr)^ir)-rEir)^Eir)^ (3) 
— oc 

The further argument proceeds in its simplest form if the Dirac Hamiltonian possesses an 
energy reflection symmetry; viz. if there exists a unitary matrix R, which anti-commutes 
with the Dirac Hamiltonian HR + RH = 0. Then R acting on negative energy states 
produces positive energy states and vice- versa, while the zero energy state is an eigenstate 
of R. 

R'^±e = '^tE, R^±e = i^^E 

R-^Q = ±-^0 ^ ' 

(The effective Dirac Hamiltonians for polyacetelyne and graphene possess this property.) 
As a consequence the charge density at energy E is an even function of E, both in the 
vacuum and soliton sectors: Tp%i^E = '4'-e'^-e, "^e^e = "^-e^ -Et^ Pe = P-E', and 
([3]) may be presented as 

oc 

p=^ I dE (5) 

— oc 

The vacuum continuum wave functions V'E iii © are complete; the solitonic vave functions 
^E are one short of being complete because the zero energy state ^'o is absent. Therefore 

p(r) = -l^S(r)^o(r) (6) 

and the charge Q is 

Q = j drp{v) = -\. (7) 

Although this simple argument does not show it, one can prove that the fraction is an 
eigenvalue without fluctuations, rather than an expectation value with fluctuations — the 
latter would not be interesting. 

One may easily contemplate Hamiltonians that do not possess energy reflection. For 
example, take a Hamiltonian with that property, and append to it a term proportional to 
the previously described matrix R. The new Hamiltonian no longer anti-commutes with 
R, and calculation of the fractional charge becomes much more involved. 

Two methods have been developed for dealing with this more difficult situation. One 
can take a field theoretic approach and calculate in perturbation theory the induced vacuum 
current in the presence of a soliton < tp^y^ip > . Here tp'j^ tp is the field theoretic current 
operator, and < | > signifies the field theoretic vacuum in an external background. This 
determines the induced charge density and therefore the charge. [The induced current 
approach is available only for problems without energy reflection symmetry — the charge 
in the symmetric case can be obtained in the limit of vanishing asymmetric effects [5].] 
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Alternatively one can show that the charge density is related to the spectral asymmetry in 
the soliton sector [B]. 



(8) 



The calculated vacuum charge in the absence of energy reflection can become an irra- 
tional quantity, approaching — ^ when the strength e of the term violating energy reflec- 
tion vanishes (But later I shall describe a more involved scenario.) Note that the spectral 
asymmetry expression ([8]) immediately yields Q = —1/2 when energy reflection symmetry 
is present, because in that case the integral vanishes. 

All these are fascinating quantum effects, but it remains a challenge to find them in 
the functional interval. 

Next I shall turn away from elaborate formalism, and give a simple counting argument 
that establishes the fractionalization effect in one dimension, in the linear polymer of 
polyacetylene. 

To begin a description of polyacetylene, we imagine a rigid array of carbon atoms, 
about lA apart, exhibiting a left-right symmetry. One might expect that thermal and 
quantal fluctuations lead to oscillations about the equidistant equilibrium positions. But 
in fact something more dramatic happens. The energetics of the system force the atoms 
to shift by about .04A to the left or to the right — both are allowed due to the left-right 
symmetry. This is a consequence of an instability of the rigid lattice, identified by Peierls. 
Therefore the polyacetylene chain presents two equivalent vacua A and B, (see Fig. [1]). 
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Figure 1: The equally spaced configuration of carbon atoms in polyacetylene (O) possesses 
a left-right symmetry, which however is energetically unstable. Rather in the ground 
states the carbon atoms shift a distance /i to the left or right, breaking the symmetry and 
producing two degenerate vacua (A, B). A soliton (S) is a defect in the alteration pattern; 
it provides a domain wall between configurations (A) and (B). 
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The potential energy of the distortion field (phonon) shows a familiar double well shape: 
the left-right symmetric point at the origin is unstable; stable configurations at the two 
minima break the left-right symmetry (see Fig. [2|). 




Figure 2: Energy density V{(j)), as a function of a constant phonon field (p. The symmetric 
stationary point, = 0, is unstable. Stable vacua are at cp = + (A) and0 = — |(/>o|, B. 

By now it is well known that a double well potential, like in Fig. [21 supports also a 
kink configuration that interpolates between the two vacua. Physically this represents a 
defect in the bonding pattern. These profiles are shown in Fig. [3l 
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Figure 3: The two constant fields, it | (/)o |, correspond to the two vacua (A and B). The 
two kink fields, ^(ps, interpolate between the vacua and represent domain walls. 
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Let us now consider a sample of the A configuration, with two sohtons — two defects — 
inserted. Let us count the kinks in the A configuration without and with two sohtons. These 
numbers need be considered only in the region between the solitons, because elsewhere the 
patterns are identical (see Fig. 




Figure 4: Two soliton state carries one fewer link relative to the no-soliton vacuum A. 

One observes that in presence of the two solitons, there is one fewer kink. Imagine now 

separating the solitons to great distance, so that each acts independently. We conclude 

that the 1-kink deficit must be equally divided, producing a fermion state with number 
_ 1 
2- 

In fact this has been indirectly observed. However, the theory has to be elaborated 
before it confronts experiment. Our argument has ignored spin. Since electrons have 
spin ±^ and two fit in each level, all our results are doubled. The charge defect in the 
presence of a soliton is —1, but there is no net spin since all other electrons are pairwise 
aligned. Inserting one electron erases the defect, but produces a spin ^ excitation with 
zero charge. This charge-spin separation, which has been observed experimentally, gives a 
physical realization to the phenomenon of charge fractionalization in polyacetylene [7] . 

The Dirac Hamiltonian matrix, relevant to this 1-dimensional problem, is a 2 x 2 matrix 
acting on a 2-component "spinor", thereby producing a matrix equation in one spatial 
dimension 

H{lp) = ap + P(p 

The kinetic term ap{a = a^,p = —idx) comes from linearizing a quadratic dispersion law 
around the two points, called "Dirac points," where it intersects the Fermi surface, thereby 
giving rise to a 2-component structure. The interaction with the phonon field fi^p (/? = o"^) 
produces a gap with homogenous tp and a zero mode with a kink profile for ip. The 
matrix anti-commutes with H{ip) and acts as the R matrix that implements the energy 
reflection symmetry, which is lost when ea^ is added to the Hamiltonian, thereby producing 
an irrational charge. 

More recently there appeared a 2-dimensional material, for which a similar analysis 
has been performed. This is graphene, which is described by a hexagonal lattice of carbon 
atoms that is presented as a superposition of two triangular sublattices see Fig. 5. Each 
sublattice supports two Dirac points, hence a 4-dimensional matrix equation on the two 
dimensional plane becomes relevant. The background field is taken to be a combination of 
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Figure 5: The hexagonal graphene lattice, taken as a superposition of two triangular 
sublattices. 

scalar and vector A [8j. The Hamiltonian governing electron motion reads 

H{ip,A) =Q;-[p-75A]+^[(/j^-i75(^^]. 

Here " = (^q , = (cr\cj2); 75 = 

where 99^, ip^ form the real and imaginary parts of a complex scalar field ip = ip"^ + i<p^ . The 
matrix that effects energy reflection is i? = = 

The topologically trivial background consist of constant and vanishing A. This gives 
rise to a gap in the energy spectrum. For the topologically non-trivial back ground we take 
for ip and A a vortex profile. A zero mode ensues, and charge becomes — ^ (with a single 
vortex). The mid-gap state is bound with just a vortex configuration for ip and vanishing 
A. A is not needed for binding the zero energy state. Its presence does affect the profile of 
the zero-energy wave function, but not the vanishing of the eigenvalue. However, a vector 
potential is needed to give the vortex finite energy: the scalar field and vector field vortex 
configurations separately carry infinite energy, which is rendered finite in the combination. 
In other words a pure ip vortex cannot "move" because it is infinitely heavy. Evidently 
this infinity is screened away by the contribution from A. 

When eR = ea^ is appended to H, the energy reflection symmetry is lost, and the 
induced charge structure becomes interestingly complicated. For a pure ip vortex, without 
a vector potential, the charge becomes irrational — as expected and as it happens in the 1- 
dimensional case. However, including a vector potential in its vortex configuration changes 
the charge back to — ^ — its value in the energy refiection symmetric case. Just as for 
energy, the vector potential screens away a contribution from the scalar field. 
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The above 2-dimensional story is very interesting in its elegant intricacies. But it must 
be stated that its vahdity as a theoretical description of actually realized physical graphene 
remains to be established. Also as yet there is no experimental verification of the charge 
fractionalization phenomenon in graphene. Nevertheless it would be most interesting to 
develop a functional integral analysis of this phenomenon. 

This work was supported in part by funds provided by the U.S. Department of Energy 
under cooperative research agreement #DE-FG02-05ER41360. 
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